Imbalanced ultracold Fermi gas in the weakly repulsive regime: Renormalization 

group approach for p-wave superfluidity 
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We theoretically study a possible new pairing mechanism for a two-dimensional population imbal- 
anced Fermi gas with short-range repulsive interactions which can be realized on the upper branch 
of a Feshbach resonance. We use a well-controlled renormalization group approach, which allows an 
unbiased study of the instabilities of imbalanced Fermi liquid without assumption of a broken sym- 
metry and gives a numerical calculation of the transition temperature from microscopic parameters. 
Our results show a leading superfluid instability in the p-wave channel for the majority species. The 
corresponding mechanism is that there are effective attractive interactions for the majority species, 
induced by the particle-hole susceptibility of the minority species, where the mismatch of the Fermi 
surfaces of the two species plays an important role. We also propose an experimental protocol for 
detecting the p-wave superfluidity and discuss the corresponding experimental signatures. 
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I. INTRODUCTION 

Much of the interest in uhracold atomic gases conies 
from their amazing tunability. Experiments on ultra- 
cold atomic gases allow fermionic pairing phenomena to 
be manipulated much more precisely and controUably 
than those in solid state systems. There are many im- 
portant experiments in ultracold gases which undoubt- 
edly illustrate this advantage, such as the crossover from 
Bose-Einstein condensation (BEC) to Bardeen-Cooper- 
SchriefFer (BCS) superfluidity with the help of Feshbach 
resonance and superfluid- Mott insulator transi- 

tions with optical lattices d, 0] . 

Due to the wide-range tunability of the effective in- 
teratomic scattering length there are strong motivations 
to study the pairing phenomena with population imbal- 
anced ultracold Fermi gases in different regimes. How- 
ever, in ultracold Fermi gases, in contrast with solid state 
systems, the pairing state is not easily achieved due to 
the smallness of the gap parameter. Therefore, previous 
investigations on superfluidity of imbalanced Fermi gases 
mostly focused on the unitary regime where the scatter- 
ing length is large [7l-[l0|. In systems with attractive in- 
teractions, the presence of population imbalance can en- 
rich the possibilities for pairing states. As predicted by 
previous works, there ma y b e Larkin-Ovchinnikov-Fulde- 
Ferreh (LOFF) state [IMI, breached pair state 
and deformed Fermi surfaces 17| . Pairing can also occur 
when there are intermediate bosons for providing effec- 
tive attractive interactions such as bosonic molecules in 
deep EEC regime, where there may be p-wave superfluid- 
ity jl8H2(| , and phonons of a dipolar condensate [2l| . Be- 
sides, in a system where the bare interactions are purely 
repulsive, there are also possibilities for effective attrac- 
tive interactions to emerge. It was first studied by Kohn 
and Luttinger 22] , where a three-dimensional (3D) elec- 
tron system was considered. In 3D electron systems, the 
particle-hole susceptibility xi^) has a strong k depen- 



dence for k ^ 2fci?, which is responsible for the emer- 
gence of effective attractive interactions in high angular- 
momentum channel. However, dimensionality can signif- 
icantly change the behavior of x(^)- In two dimension 
(2D), x(^) is momentum independent when k ^ 2fcF, 
and there may be superfiuid instability in the presence 
of population imbalance |2^ , which is different from the 
Kohn-Luttinger type. 

In this paper, we consider a population imbalanced 
2D ultracold Fermi gases in the weakly repulsive regime, 
which can be realized on the upper branch of a Fesh- 
bach resonance [l^. There are two novelties in our sys- 
tem that should be emphasized. Firstly, the bare inter- 
actions between atoms of two different hyperfine states 
are purely repulsive. Secondly, there are no interme- 
diate bosons for providing effective attractive interac- 
tions such as phonons in traditional superconductors or 
bosonic molecules at the BEC side of a Feshbach res- 
onance. Our study shows that there is an alternative 
choice of p-wave superfluid state induced by the popula- 
tion imbalance, which fundamentally breaks the spin ro- 
tation symmetry. The particle-hole susceptibility of the 
minority species can induce an attractive interaction for 
the majority species because of the population imbalance. 
This mechanism of superfluidity resembles qualitatively 
the situation in the Ai phase of superfluid ^He [1^ and 
2D electronic gases (23. [26|. 



Our theoretical framework is heavily based on the 
renormalization group (RG) theory for interacting 
fermion systems |27h30| | . The RG framework provides us 
a powerful tool to treat competing instabilities simultane- 
ously, and most importantly, to justify the leading insta- 
bility channel (28j . Furthermore, we can identify the crit- 
ical temperature from the onset of the instability channel 
(3l| . By performing RG process at finite temperature and 
solving the flow equations numerically, we can obtain the 
phase transition between normal state and p-wave super- 
fluid state. Within this framework, in the second stage 
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of RG, when mode eliminations have reached an momen- 
tum cutoff A far smaller than the Fermi momentum Kp^ 
a large-N expansion emerges with N = Kpj t^^ which is a 
strong suggestion for us to extend our results from weak 
to intermediate coupling regime .28]. 

The paper is organized as follows: The first stage of the 
RG approach for building the model of interacting imbal- 
anced fermions is described in Sec. II. Sec. Ill illustrates 
the non-perturbative RG method for unequal Fermi sur- 
faces and obtain the flow equations. The RG analysis 
indicates a leading instability in the p-wave Cooper chan- 
nel when the population imbalance is present. In Sec. IV 
we numerically solve the flow equation at flnite temper- 
ature. We obtain the critical temperature at which the 
normal Fermi liquid state becomes unstable in the p-wave 
Cooper channel. Furthermore, with the large-A'^ analy- 
sis, we extend our results from weak coupling regime to 
intermediate coupling regime where we may have higher 
critical temperature. Sec. V contains experimental dis- 
cussions and conclusions. 



II. MODEL BUILDING: THE FIRST STAGE OF 

RG 

We consider a population imbalanced Fermi gas with 
short-range Hubbard repulsive interactions, whose parti- 
tion function 

Z = y ^^]e^^[^''''l (1) 

with 

= 5o[v5,V'] + ^/[vi,V'], (2) 

where Sq is the free part, 

k is short for {kn,k). a =t or ], represents two different 
hyperflne states. E(k) = /{2m) is the free energy of 
atoms. is the chemical potential, and the Fermi mo- 
mentum satisfies /Xcr = Kp^/{2m). Population imbalance 
is put in by setting fi^ ^ fi^ [ Without loss of generality, 
we can assume fj,-^ > fi^. ). We work at finite temper- 
ature, with imaginary time formulism, where fc„ is the 
fermionic Matsubara frequency. The interacting part of 
the action reads 

= ^ X! ^t(^l)V'i(fc2)V'i(fc3)V't(^4) 

x4i„+fe2„,fc3„+fe4„'5(fci + fc2 - fca - fc4)-(4) 

The basic idea of RG is to gradually integrate out 
"faster" degrees of freedom which have larger momen- 
tums locating in a shell region in momentum space and 
see how the resulting effective Hamiltonian will flow un- 
der such process. To begin the RG process, we first in- 
troduce an artificial energy scale VIq and integrate out 



degrees of freedom with energies higher than f2o to ob- 
tain the effective theory of the system at energy scale Hq . 
We assume that the bare Hamiltonian is in the weak cou- 
pling regime. Thus, if we choose ilo not much lower than 
the ultraviolet cutoff of the bare Hamiltonian, the "inte- 
grating out" procedure can be done by using straightfor- 
ward perturbative approach because there is no signifi- 
cant renormalization of the coupling parameters. 

In detail, we first divide the degrees of freedom of the 
system into "slow modes" 

^< = Mk),^^ = Mk) for \ek,a\ < (5) 

and "fast modes" 

V'> = Mk),^b> = Mk) for \sk,a\ > f^o, (6) 

where £k,o- = E{k) - 

We then carry out the "modes elimination" by inte- 
grating out the fast modes and this can be formally writ- 
ten as 

Z = J D[Vi,^]e-^['^''''l 

^ J i^[v:;<,v<]e-^"°[*''*']. (7) 

After gathering all terms independent of ip-^ and 
into Slip^jip^] , the rest terms can be written as 
-S'[V;>,V'>] - 5*2 V'^^^], and is the result- 
ing effective action at energy scale JIq. 
Generally, has the form: 

[V;, V] = 50^" ^] + [V;, V] , (8) 

where Sq°['iP, V] is the free action, 

no 

^o""" [Vi, V-] = E ^ik)i-ik„ - ^ -t- E{k))^{k) (9) 

k 

and ij)] is the interacting action which will have the 

most generic form after the "integrating out" procedure. 

1 "° 

S^°[Tp,^]^— E u(fci,CTi,fc2, 0-2,^3, era, ^4, 0-4) 
Xl/^cri {kl)lJ>cT2 {k2)lpa3 (fc3)l/'<T4 (^4) 

^ ^ki„+k2n,k3n + k4,„ ^ C 1 +Cr2 - CTS ~ (^4, ,0 

xS{ki + k2-k3-k4). (10) 

where the superscript ilg of the summation operator 
means the summation is done within the slow modes 
characterized by the energy scale ilo (Eq. ([SJ). 






FIG. 1. Interaction vertex within the same spin species at 
order . It is induced by the susceptibility of "fast modes" 
with opposite spins. The propagators of fermions with spins 
opposite to the external ones are represented by dashed lines. 
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FIG. 2. Four-point vertex for the effective action at energy 
scale fio in the Cooper channel at one-loop order. Loop mo- 
mentums are restricted in the thin shell around Fermi surface. 



With singlet pairing suppressed by imbalance, we shall 
consider triplet pairing between fermions with the same 
spin. As shown in Fig. [T] the induced interaction within 
the same spin species in the Cooper channel is of order 
U'^ and depends on external momentums and frequencies. 
However, only the constant term of a coupling function 
is not irrelevant in the tree level scaling [2g]. Therefore, 
we will focus on the induced Cooper channel effective 
interaction of order U'^ with external momentums set on 
the Fermi surface and external frequencies set to zero. In 
this case, the interaction vertex will only depend on the 
orientations of the incoming and outgoing momentums 
and can be written as 



T,{k,k) = U\-,{k~k'). 
Xo-(fc) is the susceptibility at zero frequency: 

X<y{k)= I G„{ipn,p)G„{ipri,P + k) 
Jp 



(11) 



^p-\-h.(T ^P,(7 



(12) 



where / is the Fermi-Dirac distribution function, and 
Ip = Sp„ / (tJ)'^ ■ Because the susceptibility is not 
singular in the limit JIq — >■ (Oq ^ /i), we can only keep 
the zeroth order term by setting — 0. 



III. THE SECOND STAGE OF RG 

After obtaining the effective action at energy scale fto 
around the Fermi surface, Shankar's RG [l^ for fermions 
can be carried out. We use the field-theory approach 
and calculate the four-point vertex at one-loop order at 
energy scale V, (see Fig. [5]): 

rW(fe,fc') = r.(fc,fc') 

T„{k,p)V„{p,k') 



{iPn 



r)(-«Pn 



(13) 



where means the momentum integral is restricted in 
the shell region in momentum space with energy devia- 
tion less than f2 with respect to Fermi energy, and f2 is 
a energy scale within fio, i.e. < fio- The one- loop 



correction can be further carried out as: 



■E 



pAp 



2-K 



2-rr 



de T,{kp)T,{p,k') 

27r [ipn - ep,a){-ipn ~ £p,a) 



■ pdp 
"27 



2TT 



de r,(fc,p)r,(p, k') ^^^^^^ ^ 



27r 



2e, 



(14) 



where is the momentum cutoff corresponding to fl, 
and / " is short for /|p_x^^|<A„- to 17 < fi^), 

we can approximate Q as vpa-Aa, where vpa is the Fermi 
velocity of atoms with spin a. 

Because the four-point vertex is related to the scatter- 
ing amplitude of certain scattering process, which is a 
physical observable, it should not depend on cutoff: 



dr^"^(fc,fcO 
di 



= 0, 



(15) 



where I = ln(17o/f^) = ln(Ao/A). 

From Eq. ([T5|) . we can get the flow equation as follow- 



ing 



dr^(fc,fcO 

dl 



-Pj — r^(fc,p)r„(p,fe )tanh — 



(16) 

where p — m/2n is the 2D density of states. This is an 

ordinary differential equation for matrix Ta{k,k') with 
initial condition: 



r,(fc,fc';17o) = C/2x-.(fc-fe'), 



(17) 



where Ta{k, k' ]VIq) denotes the effective interaction ver- 
tex at energy scale fio- For convenience, we de- 
fine a dimensionless coupling function g^ik^k' = 
pTa{k,k';D.o). In the presence of rotational symmetry, 

gcrik, k') only depends on the relative angle between the 
incoming and outgoing momentum, and the flow equa- 
tion can be decomposed into uncoupled equations for 
eigenvalues of channels with different angular momen- 
tums [301: 



dA„ 



= - {Km) tanh— , 



dl ' 2 

where m labels different angular momentum channels 



(18) 
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FIG. 3. Induced interaction vertex at energy scale fio after 
reconsidering a order correction, which is responsible for 
removing the flo dependence. The loop modes are the "fast" 
ones, whose corresponding on-shell energies are larger than 



The right hand side of Eq. (|T8|) is negative definite, 
which means that for an initially attractive channel, 
Xam may be renormalized to negative infinity as the en- 
ergy scale goes down to the Fermi energy. A qualita- 
tive argument of the critical temperature can be given 
based on Eq. (|T8)) . At low temperatures, tanh(/3r2/2) 
equals almost unity for nearly all fi's when > 0, and 
drops rapidly to zero as fl approaches zero from about 
n ksT. Therefore, when temperature is low enough, 
we can approximate tanh(/3f2/2) to be unity, and easily 
get the solution of Eq. ([TOl) . which guarantees a diver- 
gence at a certain energy scale. If this energy scale is 
higher than ksT, the approximation used above is self- 
consistent, and the divergence indicates the superfluid 
instability. In other words, this energy scale gives a qual- 
itative estimation of the critical temperature. 

Results obtained by using the methods introduced 
above will certainly depend on the energy scale JIq. How- 
ever, it is clear that J7o is more a calculation device than 
a physical energy scale [1^. Any physical predictions 
should not depend on flo- In fact, like RG in quantum 
field theory, we can make results independent of fJo a-t 
any order of U. This can be achieved by simply consid- 
ering diagrams of the same order as the ones in the second 
stage when carrying out the perturbative calculations in 
the first stage [2^, Is^] ■ 

In detail, we have to take an additional term into ac- 
count (see Fig. [3]) and take Tc^ik, k' ;ilo) as the initial 
condition for the second stage. 



f,(fc,fc';rio) 



which has a similar form as Eq. except that the re- 
gion of the momentum integral is the area besides the 
thin shell around the Fermi surface and the dependence 
of the vertex on the magnitudes of momentums and fre- 
quencies should also be considered. 

We first consider the frequency summation in the sec- 
ond term of Eq. Most contribution comes from 
neighborhood of po ~ and \p\ ~ Kpa-, so we can first 
set frequencies in T's to zero. Then the frequency sum- 



mation can be carried out as: 
f,(fc,fe';r!o) 



r.(fc,fc')- 

X tanh ^£pf^ . 



pdp deT^{k,p)T,{p,k') 



2tt ./(, 2tt 



2e 



per 



(20) 



where A^-o = ^o/'i^Fa- 

Because the bare interaction is short-range and the 
Fermi gas is dilute which means the interatomic distance 
is much larger than the range of the interatomic inter- 
action i.e. Kp^ ^ ro, we can introduce an ultraviolet 
cutoff Ah = 2Kf\ here. In the second term of Eq. (pn)) 

does not vary much with respect to momentum within 
this cutoff f23j at low temperatures compared with the 
notable dependence of the other term on momentum. 
Therefore, we can neglect the dependence of the F's on 
the magnitude of the momentums. Using the dimension- 
less coupling function which has been defined above as 

gcr(fc, fc'; i7o) = pVa{k, k'-jfto) the initial condition of the 
fiow equation can be written as 



5^(fc,fc';fio) 
9aik,k') - 



d0 



g„{k,p)g„{p, k') 



/o 2^ 

where F„(yi) is an auxiliary function defined as 



pjk 



pdp 1 



■ tanh 



I3e 



pa 



(21) 



(22) 



With the help of Eq. 
pressed as: 



ep„|<n 27r 2epa 2 

, fiow equation can be ex- 



dg„{k,k') 
dil 



de 



^ ^^9<y{k,p)g„{p,k)- 



(23) 



It can also be written in a more compact form, regard- 
ing g„{k, k') as a matrix {ga)},^^, 



dn 



dF^jn) 

dn ' 



(24) 



As we mentioned earlier, because of the rotational sym- 
metry of the system, the coupling function can be de- 
coupled in the angular momentum representation, and 
Eq. (p4)) becomes a series of flow equations of individual 
eigenvalues, 



dA„ 



dfl 



.dF, 

dn 



(25) 



with innitial conditions, 

Xami^o) = Xam ' (A,™)^ (F^inn) - F,(l]o)) , (26) 
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which can be easily integrated and gives 

-F{n)+Fino). (27) 

To order [/^, we have 

X^Un) = (A-^ + F^inn) - F,{n))-\ (28) 

which is independent of flo, and controls the flow be- 
havior of the coupling strengths. At zero temperature, 
as mentioned above, a negative eigenvalue will flow to 
infinity at certain energy scale and cause superfluid in- 
stability. As temperature goes higher, the divergence will 
appear at lower energy scale. When temperature reaches 
a certain critical value, the divergence will not arise until 
we renormalize to the Fermi surface. Above this critical 
value, no divergence exists during the whole renormal- 
ization process down to the Fermi surface, which means 
that the Fermi liquid state is stable in the corresponding 
channel. 



IV. NUMERICAL SOLUTIONS OF RG FLOW 
EQUATIONS 

In this section, we will determine the critical temper- 
ature by solving the flow equations numerically. As ex- 
plained above, at critical temperature, we have 



= F,(17ff)-F,(0). 



(29) 



For convenience, we can define some dimensionless pa- 
rameters as follows Xa = k/ Kp„,tc! = TjTp^^ where Tp 
is the Fermi temperature {Tp — iialks), and rewrite the 
susceptibility and the eigenvalues to the final form for 
numerical calculation. 

Susceptibility only depends on the magnitude of mo- 
mentum and can be written as 

(27r)^ Jo Jo ^p. 



y 



1 



(30) 



Eigenvalues of the dimensionless coupling function ga 
have the form 

Xarn^pU'^ — X-ct (2fcF£T siu - ) C0s(to6') . (31) 

Jo 27r 2 
For the p-wave and /-wave case, we have respectively 

y 1 - 2a;2 



- -{pUf- f dx [ dy^ 
I" Jo Jo \/l 



1 



(32) 
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FIG. 4. (Color online) Strength of the pairing interactions ( 
—Act/ (pU)^) for spin spieces a in the p-wave (top) and /-wave 
(bottom) channels at different temperatures, where rjcr is the 
imbalance ratio defined as rja = Kf-<j / Kpa- The scales of 
the vertical axes of these two figures are chosen to be the same 
for convenient comparison. 



\a^ = -{pU?- ['dx f dy JL^ 

TT Jo Jo VI - y 

-32x^ + 48x4 - 18a;2 -f 1 



x^ 



1 



(33) 



First, the eigenvalues of p-wave (Acri) and /-wave 
(Ao-s) channel are obtained at different imbalance ratios 
[rjfj = Kp-a / Kpa) and temperatures as illustrated in 
Fig. m where it can be seen that eigenvalues of the p- 
wave channel is more negative than those in f-wave chan- 
nel, indicating the leading instability in p-wave channel. 
Eigenvalues are almost zero when rj^ > 1, which means 
that there is no obvious instability for smaller Fermi sur- 
face (we will focus on the majority speices in the remain- 
ing part of this paper). This is similar to the Ai phase of 
■^He [1^ when applied a magnetic field which will cause 
spin population imbalance. According to the qualita- 
tive arguments on phase given by Leggett [33[, the 
reason why pairing only happens for the bigger Fermi 
surface is that the density of states at the bigger Fermi 



0.4 0.5 0.6 0.7 0.8 0.9 1.0 
pU 

FIG. 5. Critical temperature (Tc) from weak coupling regime 
to intermediate coupling regime with imbalanced ratio (ry) set 
to 0.5 near the optimal value. Tp is the Fermi temperature. 
The variation of Tc with respect to rj at different coupling 
strengths will be illustrated in Fig. (6] 

surface is larger, which results in higher critical temper- 
ature for the majority species. However, in a 2D system, 
density of states is a constant and we are thus facing a 
different situation from the Ai phase of '^He. Besides, 
we can see that, for p-wave channel, the optimal imbal- 
ance ratio where the most negative eigenvalue appears is 
about 0.5, which is consistent with Ref. [1^. When tem- 
perature goes higher and becomes comparable with the 
Fermi temperature, the eigenvalues are concealed under 
thermal fluctuation. 

To guarantee the validation of the perturbative ap- 
proach in the first stage, the dimensionless coupling pU 
should be small. The induced vertex, which is of order 
[/^, will be even smaller, and exponentially suppress the 
critical temperature. Considering the large-N emerging 
in the second stage of RG [2I], we can extrapolate the 
results to the intermediate coupling regime [23i|. 

Setting imbalance ratio to 0.5, near the optimal value, 
we plotted the critical temperature from weak to inter- 
mediate coupling regime (see Fig.[S]). It can be seen that 
the critical temperature drops quickly as the interaction 
strength goes smaller. 

We also calculated critical temperatures at different 
imbalance ratios with coupling strengths pU ^ I in the 
intermediate regime (see Fig. [5]). For fixed coupling, the 
highest critical temperature appears near 77 = 0.5 as ex- 
pected and is about lO^^Tp but also drops quickly as the 
imbalance ratio tends to zero or unity. 

V. DISCUSSIONS AND SUMMARIES 

We have shown that the population imbalance induced 
p-wave superfluid state may be observable in a 2D repul- 
sive fermion gas. Population imbalance can be achieved 
by an unequal mixing of atoms in two hyperfine states, 
and tunable repulsive interactions can be realized by 
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FIG. 6. (Color online) Critical temperatures Tc in unit of Tf 
at different imbalance ratios (rj) and coupling strengths (pU) 
in the intermediate coupling regime. The variation of Tc with 
respect to pU from the weak to the intermediate coupling 
regime is shown in Fig. [5] 



usin g t he upper branch of a Feshbach resonance. In 
Ref. |24| ■ ^Li atoms in the repulsive regime were used to 
study the itinerant ferromagnetism. One problem that 
should be considered is that the upper branch of a Fesh- 
bach resonance is an excited branch, and will decay to the 
BEC molecule state due to inelastic three-body collisions 
[33 |. However, with small scattering length and popula- 
tion imbalance, the decay rate is suppressed [13, and 
the system may be metastable for observation. For ex- 
perimental observations, we suggest to look for rotational 
asymmetries in the momentum distribution or pairwise 
correlation in the time of flight expansion images of the 
dominant species [s^, Is^- In addition, the transition 
temperature can be raised in two different ways (see 
Fig. [5]). One is to adjust the imbalance ratio to the opti- 
mal value which is around 0.5. As can be seen from Fig.[5l 
the theoretical transition temperature should be around 
lO^^Tp ^ lO^^Tp in the weak repulsive regime near the 
optimal imbalance ratio. Another one is to increase the 
coupling strength by Feshbach resonance. By extending 
our result to the intermediate coupling regime, we get 
an estimation of the critical temperature which reaches 
as high as 10~^Tp. Since our approach is asymptotically 
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exact, the perturbative calculations are well controlled in 
the first stage of RG. After safely arriving at the second 
stage, where the cutoff A is much smaller than the Fermi 
momentum Kp, the emergence of a large-N ensures the 
non-perturbative nature of the momentum shell RG in 
the second stage where a quantitatively calculation of 
critical temperature was given. However, for comparing 
with the experimental results, we should notice that an- 
other important issue is the trap effects on our system. 
In striking contrast with s-wave superfluid state, the trap 
asymmetries would have a strong influence on the sponta- 
neously preferred orientation of p-wave superfluid state. 
We will study the trap effects in our future work. 

In summary, we studied a possible new superfluid state 
for a 2D population imbalanced fermion gas with short- 
range repulsive interactions. This phenomena is different 
from the ones in the BEC-BCS crossover where the BEG 
molecule state is concerned. It is also different from the 
ones in the unitary region where the scattering length 
is approaching infinity and many universal properties 
emerge. For the system considered in this paper, the bare 
interaction is purely repulsive, and there are no inter- 
mediate bosons for inducing attractive interactions. We 
studied this system based on RG approach and treated 
different instabilities on an equal footing. There are no 



assumptions of specific orders compared with mean-field 
approach. What is essential for the Cooper instability 
is the mismatch of the Fermi surfaces caused by popu- 
lation imbalance in our system. It can also be achieved 
in mixture of fermions with unequal mass, to which our 
approach can be generalized straightforwardly. By work- 
ing in the finite temperature formulism and numerically 
solving the fiow equation, we gave a quantitatively calcu- 
lation of the critical temperature. Our study is of partic- 
ular significance both for probing p-wave superfiuidity in 
the novel regime experimentally and for studying imbal- 
anced fcrmionic systems with RG theory theoretically. 
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